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(Abstract.*) 

Introduction. 

The mean-free-path phenomena of gases can be explained in a general 
way by a very elementary form of the kinetic theory, but to develop a 
satisfactory mathematical treatment .of them, yielding accurate numerical 
results, is less easy. The difficulty varies to some extent with the nature of 
the molecular model chosen as basis. It is curious that the adoption of one 
particular model, viz., a centre of force varying inversely as the fifth power 
of the distance, removes nearly all the analytical complications in the theory. 
The first accurate treatment of viscosity, conduction, and diffusion, due to 
Maxwell, was made possible by this fact, and it is remarkable that Maxwell's 
theory requires no knowledge of the velocity-distribution function corre- 
sponding to the non-uniform state of the gas. If the molecules are of any 

* This abstract summarises two papers in the 'Philosophical Transactions/ the first 
(received October 5, 1915) dealing with simple gases (' Phil. Trans.,' A, vol. 216, 
pp. 279-348, 1915), and the second (received May 19, 1916) with composite gases. These 
two papers amplify and complete an earlier memoir (' Phil. Trans.,' A, vol. 211, pp. 433- 
483 (1911)), which contains a first approximation to the present theory. 
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2 Dr. S. Chapman. The Kinetic Theory of 

other type, the function must be determined, a task wherein lies the main 
difficulty of the investigation. It is, however, of more than mathematical 
Interest to widen the basis of the theory, since Maxwell's molecular model 
does not satisfactorily represent the molecules of actual gases. 

The velocity-distribution function has been proved by Boltzmann to 
satisfy a certain integral equation, which Hilbert has shown to be sufficient 
to determine the function uniquely. Lorentz has solved an important 
dege aerate form of this equation, while Pidduck* has obtained solutions of 
some typical examples of the general case, and so deduced many interesting 
numerical results relating to diffusion. The present method, however, is not 
based on Boltzmann's equation. Though less direct, it seems to possess 
advantages not shared by the other method, in that it yields results of 
complete formal generality, which can be reduced to numerical form with 
comparatively little arithmetical labour. To comply with the requirements 
of the pure mathematician, the results should be proved to be in conformity 
with Boltzmann's equation. Progress has already been made in this 
direction, and has brought to light several theorems of much analytical 
interest. As the insertion of these would be unsuitable in memoirs intended 
primarily to deal with questions of physics, they are reserved for a future 
paper. 

Part I. — Outline of the Mathematical Method. 
(1) Notation and Analysis of the State of the Gas. 

We consider a gas composed of a mixture of two sets of spherically 
symmetrical molecules of masses mi, m 2 ,f in numbers v\, v 2 per unit volume 
at (x, y, z, t).l The component densities p X) p 2) are consequently equal to 
vi'irii, v 2 m 2 . The external force acting on each molecule oii\ (and similarly 
for the molecules m 2 , with change of suffix throughout) will be denoted by 
(Xi, Yi, Zi), or, vectorially, by P, while the mean velocity of this group of 
molecules at (x, y, z, t) will be denoted by (u h v\, w{) or vectorially by c\. The 
proportions of either kind of gas will be denoted by Xi, \ 2 , so that 

where vo ~ v\ + v 2 ; evidently \i -f \ 2 = .1. 

* Boltzmann, ' Vorlesuugen tiber Gastheorie,' vol. 1, p. 114; Hilbert, 'Math. Ann, 
(1912); Lorentz, 'Theory of Electrons,' p. 268; Pidduck, 'Proc. Lond. Math. Soc, 
(2), vol. 15, p. 89 (1915). Also cf. Lunn, 'Bull. Amer. Math. Soc./ vol. 19, p. 455 (1913) ; 
finskog, ' Phys. Zeit.,' vol. 12, p. 58 (1911). 

t For convenience we shall suppose the heavier gas to be taken as gas (1), so that 
,w,j > m 2 . 

J I.e., at the point (#, y, z), and at time t. 
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The suffix here, and throughout, refers to the composite gas as a whole. 
Thus the mean mass m%, or external force P (or X , Y , Z ) per molecule, 
and the mean velocity c (or u% t v , w ) of the gas as a whole, are given by 
equations of the type 

m® = Xim x -f \ 2 m 2 , X = X1X.1 -f A-2X2, 

and so on. Evidently, the total density p 0> equal to p\ -f p2, is also given by v^m^ 
The suffix 0., together with a dash ('), will relate to certain multiples of 

the difference between the values of various data for the component gases. 
Thus we define X ', m$\ p&, P/, c { / by the equations 

2V = /W—Xs, W = X1X3 (mi -- m 2 ), p ' = VoWIq' /\{K2 = ^o(%- w 2 ), 
P ' = XiX 2 (m 2 Pi/m — mi T 2 /m ), £ ' = XiX 2 (ci — c 3 ). 

By inverting the above equations, we obtain the following expressions for 
the original quantities m 1} Pi, etc., in terms of m , mo', P , and so on: — 

v% = mo + '/y^Z/Xi, Pi = miVo/mo + Vo'/Xi, C\ = Co + Co'/\i, 
m 2 = m^ — mo'/Xi, P 2 = m 2 P /m — Pq' /X 2 , c 2 = c () ~~Co/X 2 , 

It may readily be seen that the motion of the gas can be analysed into 
«(«) a steady motion of the gas as a whole with velocity e , together with (&) 
a motion of interdiffusion in which equal numbers (i/oft/ P er second) of the 
two groups of molecules are transported in opposite directions, the mean 
velocities of the two streams being c</f\i and — cq j\% The momentum of 
the common motion (a) is p & cq, while the resultant momentum of the motion 
of interdiffusion (due to the difference of molecular mass) is po'c ' ; it is 
•easily seen that, with the above definitions, 

piCi H- p 2 c 2 = po^o + po'co'. 

As regards the resolution of the forces into multiples of P and P ', the 
first terms, miP /w?e on mi and w 2 Po/m on m 2 , represent forces which will 
impart a common acceleration Po/mo to each group of molecules (we may 
.suppose this to modify the common stream-velocity c ). The remaining 
components of Pi and P 5 , when summed up over the v\ or v 2 molecules of 
the corresponding group, afford equal and opposite total forces ± *>oPo'; 
such interdiffusing groups of molecules as we consider in (b) will exert 
equal and opposite forces on one another, and the forces acting, on the two 
groups must be equal and opposite if the motion of interdiffusion is to be 
maintained, or modified without imparting any common velocity to them. 
As for XoYit is introduced for the sake of symmetry: D denoting any 
differential operator, clearly 

DXi = -DX 3 = p)(Xi~X 2 ) = DX ( /. 

B 2 



4 Dr. S. Chapman, The Kinetic Theory of 

We next consider the motions of individual molecules. The velocity of 
a typical molecule mi will be denoted by (tt) h (v)i, (u% or (c)i, when referred 
to the co-ordinate axes, or when referred to axes moving with the velocity 
c appropriate to the given values of (x, y t z, t), by Ui, V x , Wi or Ci. Thus, 

Ci = (c)i—cb, C 2 = (c) 2 — c , 

the second equation referring to the molecules m 2 . 

The mean value of any function of the molecular velocities will be 
denoted by a bar placed above the expression for the function. Hence, by 
definition, 



(C)i = Co + Ci = c lt (c) 2 = c + C 2 = c 2 , 

or 

Ci = Ci— c = c '/\i, G 2 = c 2 —c == — c 7^2. 

So far, C has denoted a vector quantity, and the equations involving it have 
been vectorial : henceforward it will denote the amplitude of the vector, so* 
that C is essentially positive, and 

d 2 = IV 4- Vi 2 + Wi 2 , C 2 2 = U 2 2 + V 2 2 + W 2 2 . 

The mean energy of the motion (U, V, W) " peculiar " to each molecule 
is clearly -|mC 2 per molecule, and we write 



mid 2 = 3/2Ai = 3ETi, w 2 C 2 2 = 3/2A a = 3BT 2 , 

where R is the universal gas-constant. These equations define A and T: 
Ti, T 2 , will be termed the " temperatures " of the component gases,, which 
are not necessarily exactly equal, though all the departures of the gas from 
the uniform state will be supposed small. The mean hydrostatic pressures 
of the component and composite gases, and also A and T (the " temperature" 
of the composite gas) are defined by 



p Y = ^miCi 2 = i/i/2Ai = K1/1T1, p 2 = ±v 2 m 2 C 2 2 = v 2 j2h 2 = Rz> 2 T 2 , 

jpo == ^1+^2 = vi/2hi + v 2 /2h 2 = vof2h (} = Bi/ T . 

Clearly . T = XiTi -f X 2 T 2 , 1/ A = \i/ Ai + X 2 /A 2 . 

The six pressure components p xx> p yy , p zz , p^j, p yZy p zx > for the component and 
composite gases are similarly defined, e.g., 



(?**)<> = (Pxx)l + (j>xz)i = W^l U \ + ^% U 2* 



(^y)o = (p^)i + (^)2 = ^imiUiVi4-z/ 2 m 2 U 2 V 2 . 
We also define To', po', A ', by the equations 

To' = X 3 X 2 (Ti-T 2 ), 1/V = XiX 2 (l/Ai»l/A 2 > 
po' = v /2A ' = XiX 2 (^i/Xi— _p 2 /X 2 ) = R% r iV. 
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Evidently T x = T + To'/V, T 2 = T - TV/** 

and To' is a measure of the difference between the temperatures of the com- 
ponent gases. 

(2) The Velocity-Distribution Function. 

When the gas is uniform, so that c 0) X, h are constant, while c^, h \ P, are 

zero, the function /(U, V, W) representing the distribution of the components 

of velocity (TJ, V, W) among the molecules m x or m 2 assumes Maxwell's well 

known form : 

[/i(U 1? V b Wi)] = (horK/iry&e-W>*, 

[/s(U a , V a , W 2 )] = {hm 2 jirflH^^\ 

These clearly satisfy the necessary conditions 

f [[/OX V, W)<ZUtfVdW = 1, [ [[/(IT, V, W)C 2 dUdVdW = 3/2Am. 

For the velocity distribution function in the slightly disturbed state we 
shall write 

A (U,, Vi, wo = [/!]« [l-^CUj, Vl w,)], 

where <£i is a function, of the first order of smallness, which remains to be 
determined. 

(3) The Equations of Transfer. 

Let Qi be any function of the velocity components (tt) ly (v) h (w)u of a 
molecule mi ; then z/iQi is the aggregate value of Qi summed over all the 
molecules in unit volume. The equation expressing the analysis of the rate 
of change of *>iQi is called the "equation of transfer" of Qi. It may readily 
be shown to be* 



Ct x , ij, z 



v\ v /oQi 






The rate of change is here analysed into (a) the part AQ X due to the 
encounters between the molecules m x among themselves or with the others 
7n 2i (b) the part due to the passage of molecules into or out of the volume 
considered, and (3) that caused by the action of external forces. 

By assigning to Q the values unity, (u), or (c) 2 , we may obtain the equations 
of continuity, momentum, or energy, for the composite gas, since in these 
three cases, by virtue of the conservation of mass, momentum, and energy, 
we may eliminate AQ. 

Three special functions Q are important in the present theory. By means 
of the three equations just mentioned, and by the neglect of all second-order 

* Of. Jeans' ' Dynamical Theory of Gases,' § 335 (2nd ed.). 
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quantities, the corresponding equations of transfer may be reduced to the 
following form :• — 

i.3...(2s+3) miAUlU ~»*& +*>»*&> 



m 2 AU2C2 2 * = V-vo^o' +Es^2 



1.3...(2s + 3) * *&?' 

(2ft mi) > AC o^ = cV = _ (2fe (> m 2 f A p ^ 
1.3...(2s+l) * ° 3^ 1.3...(2* + 1)' 2 

3(2& m 1 y +1 ATT2H2, _ v 3V, 2 , 9 r 
3(2a m 2 ) s4 ' 1 ATT2Q2. _ _, /A ^L4_ 2 f2sM~5W>r 

3(27i mi) g * 1 1 AywH^ _ 2 9 k^ __ 3(2/^m 2 ) g+1 1_ * y w n 2* 
1.3...(2s + 3)* lWlU ■~15 ( ^ + 0) ^™* 1.3... (2, + 3)^2 22i ' 

The following equations define f ', ^ %, and there are, of course, similar 
quantities 770', §/, %, and so on, which are defined in like manner : 

2Ao Cm? W% Cte 

c^g = 2d'ii /d%~~~dvoldy~dw /dz, c yz =s f (9%/c2? + <W/dy). 

In each of the above equations of transfer, as written, a multiple of AQ is 
equated to certain " external " or mean data of the gas. The reduction to 
this form requires no knowledge of the small deviation from Maxwell's 
distribution law of molecular velocities, i.e., of the function <£(IT, V, W), but 
such knowledge is necessary for the calculation of AQ. Conversely, however, 
the above equations for AQ throw light on the form of <£ (IT, V, W) ; in fact, 
if we assume that (j> can be expressed as a power series in IT, V, W, we may 
in this way* infer that it will have the following form : — 

4>i (Ui, YiWO == 1 2AoWaAo (Uif ' + VrV + WigiO 1\ (C1 3 ) 

+4 2/iom^o (Ui3To/ai!+ ViBTo/By + WiBTo/a^) Gi ((V) 

+ A 2A,m 1 C (aii 2 + . . • + 2c 1JZ Y 1 W l + . . . ) ■ H x (d 2 ) + D 3\o7^ • J i (Ci 2 X 

The function <f> 2 is obtainable by changing the suffix 1 into 2 throughout. 
We express the functions F(C 2 ), etc., in power series which are conveniently 
written as follows : — 

* The details of the deduction will be found in the original papers, where it is shown 
how AQ depends on <£ (IT, V, W). 



Simple and Composite Monatomic Gases. 7 

* (Cl } ~ ,r o a " 1.3...(2r + 3) Cl ' Gl (Cl } ~ ,?o A "l.3...(2r + 3)r Cl ' 

Hl (Cl } ~ Jo 7 " 1.3...(2r + 5) Cl ; J > (°i > - 2 ^ i73~(2^TI) Cl ' 

The dash ( ' ) after the sign of summation in G (C 2 ) is to indicate that r is 
to be. omitted from the numerical factor of the first term (r = 0). These 
factors are arbitrary and not really necessary, since the coefficients a, /J, y, S 
are unknown, but their insertion is analytically advantageous. Similarly it 
is convenient to denote the coefficients in F 2 (C2 2 ), etc. (which expressions 
correspond to Fi(Ci 2 ), etc., apart . from the natural change of suffix) by 
a_ r , /3_ r , 7- r , S_ r ; it is then necessary to distinguish between a r and «- n 
and so on, even when r = 0. 

The factors A , B , , D , in <£ can be chosen at will, after which the 
coefficients a, /3, 7, S, become definite ; these remain to be determined. 
Certain relations between them, however, hold good, in virtue of the 
conditions 

|Jj7(U, Y, y/)d\JdVdW = 1, JJJ/(U, V, W)C 2 d\J dYdW = 3/2Am, 

|||/(U, V, W)XJ-dUdVdW = w-«„, 

which (with the proper suffixes appended) /i(Ui, Vi, Wi) and /^(Ua, V 2 , W2) 
must clearly satisfy. These are readily seen to yield the equations : 



CO 00 



«__> ©*■ : — — v O — ?• — *• v . 



■r=0 r=0 

00 CO 



Xi 2 (2r + 3)S r = -X 2 2 (2r + 3)S-., = -2,T n '/I>,>(dX«'/dt), 

r=0 >" = 



00 00 



iXxfAof,' 2 «,.+ Bo(3T /a«) 2' r- 1 A-} 



x> 00 



In the last equation, £</ and dT /dx are entirely independent quantities ; 
hence, if we eliminate uo', we can separately equate their co-factors to zero, in 
the resulting equation. In this way the above equations may be re-written 



00 


00 


Xl S «r = 


— X2 S «— r EE 









QO CO 



Xi 2' r" 1 A- = -X 2 2' r" 1 /?-,• = - #,', 




00 00 



2Xi S rS r = —2X2 S ?'8_ r rz —&/, 
1 1 
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which also define a</> fto', and So' ; consequently, also 

* ' = i. (* ' A £ ' + /So' Bo 3T /3^), 
To' = ^8q DoTodXo /St. 

(4) Determination of the Coefficients a, /3, 7, 3. 

By means of the above expression for /(U, V, W) it is possible to 
calculate A Q, the operation (a lengthy and elaborate one) consisting in the 
evaluation of an octuple integral involving the eight variables which are 
required to specify an encounter between two molecules. Small quantities 
of order higher than the first are neglected, and the resulting expression 
for AQ always consists of a linear series of the coefficients ot r ,/3 r ,yr>Sr, 
multiplied respectively into £ '> 3T /3«r, c xx , d\o'/dt, etc. ; the factor of each 
term in these linear series depends only on v, h Q , and m, and involves the 
law of molecular interaction during encounter. It is found possible to 
obtain perfectly general expressions for these factors, i.e., no particular law 
of interaction has to be assumed. 

On inserting the values of AQ, so calculated, in the above equations of 
transfer, we may separately equate the parts on the two sides which contain 
^o', dTo/dx, c xx , or d\o'fdt. We thus obtain a number of equations, each of 
which is linear in one set only of the coefficients a r , /3 r , <y r , °r $ r > the factors 
being as above stated, viz., definitely known functions of v, ho, and m. 
There is an infinite sequence of such equations for each set of coefficients, 
one equation corresponding to each value of s, which takes all integral values 
from to go . If we solve them in the same way as for p linear equations 
in p unknowns, although p is here infinite, we arrive at solutions typified 
by the equation 



<X r = 



V \&mn) 



here a mn denotes the coefficient of a m in the nth equation of the sequence of 
equations for the «'s (i.e., s = n), and y?(a mn ) denotes the infinite deter- 
minant whose general element is a mn , while TJ r (ci mn ) denotes the same 
determinant with the elements of column m = r all replaced by unity. 
These determinants are infinite in both directions, covering the whole plane ; 
in each of the equations whence they are derived, m ranges from -co to 
-f 00 , while the equations of transfer for AQ 2 provide a second infinite 
sequence of equations corresponding to the range of n from — to — 00 . If 
we are considering only a simple gas, so that v 2 ~0, $2 and AQ 2 do not 
appear, and the determinants reduce to the more ordinary form, covering 
only a quadrant of the plane. 
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These general formulae for the coefficients a, /3, y, S, complete the determi- 
nation of the velocity distribution function / (U, V, W). The mean value 
of any function of U, V, W, can hence be calculated in the form of a linear 
series of these coefficients. It remains only to add that any such linear 
series in the «'s, 13% y% or S's can readily be transformed into the quotient 
of an appropriate single determinant by v (a mn ), v (b mn \ V (c mn ), or v (4») 
respectively. 

Part II — Application to the Theoky of the Mean-Free-Path 

Phenomena. 

(1) Viscosity. 

By I (1) and I (2),the total pressure components (p xx )o, (pxy)o> are given by 

4 1 °° 

(Pxx)o = FiWliUi 2 + V 2 m 2 U 2 2 =P0— -ftp 7TT C <> { S (vijr + V 2 J- r )}Cxx, 

4 1 



(Pwh = v 1 m 1 V 1 V 1 + y 2 ?n 2 U 2 \ T 2 = — g=g 2A"C {2(vi7 r + V27-r)}<V 

Keinembering the significance of c xx and c :qf) and comparing these equations 
with the equations of pressure of a gas whose coefficient of viscosity is k 12 , 
viz., with 

( Pxx)o =P0~ | * 12 Cxx, ( p X!/ ) = — 1 K12 C xy , 

it is clear that the composite gas behaves like a viscous fluid, and that 



2 Co 



co 



Ki2 = nn7 HT~ 2 (Vi7r + I'SY-r). 

(2) Diffusion. 
From I (3), remembering the significance of (•<>', it appears that 

— ij* Ao«o (oT"~5 ^o — tt- } + -$ *>oPo -*r~* 

\2/i da? . v m da? / * d^ 

Hence (taking the right-hand terms in order) diffusion is produced by (1) a 
concentration gradient, or variation in the relative proportions of the 
constituent gases ; (2) by external forces acting unequally per unit mass 
on the two sets of molecules, and by variations in (3) the total pressure, 
or (4) temperature of the composite gas. If we write 

D12 = — -J-AolMW, Di2' = V 12 /p 0> Dp = (mo'/mo)Di2> E>t = —^BoPo'To, 

then 

uo == —V 12 d\ ' /dx+Di 2 'voX ' + & p .dp /p dx—D T .dT /T dx, 
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and we may call Di 2 , Di 2 ', D^, and D T respectively the coefficients of diffusion,, 
forced diffusion, pressure diffusion, and thermal diffusion. The definition of 
D12 agrees with that usually given for the coefficient of diffusion. The other 
coefficients seem to be defined here for the first time. 

If, as we have supposed, the molecules m-i are the heavier, the four coeffi- 
cients of diffusion are positive. Hence, in the case of pressure diffusion, 
the heavier gas will tend towards the direction of increasing pressure ; in the 
case of thermal diffusion, the heavier gas will tend towards the direction 
of decreasing temperature. Evidently both pressure and thermal diffusion 
will in general take place during the passage of sound through a composite gas 
such as air, the effects presumably being comparable with those of viscosity 
and conduction which have . been investigated by Stokes, Kirchhoff, and 
Kayleigh. 

(3) The Transfer of Energy by Conduction, Diffusion, Etc. 

The equation of energy is the equation of transfer of -JmC 2 . On 
combining the separate equations for the two gas-components, the equation 
of energy for the whole gas, correct to the first order, is found to be 

. n DTo 9 n T (dtiQ dv (W 
Dt \ ox ay cz 



j-*~'{(*-*SM x »— w 

J L \ ox I \ ay vz / ' \ ox 



1 v 3 



Here G v is the specific heat at constant volume, and J is Joule's 
mechanical equivalent of heat. The right-hand terms, it will appear, depend 
on viscosity, diffusion, and conduction. If we neglect these small effects, the 
equation may be written 

1 DTp , o (duo , 9flo , jW\ __ 

To'd^ T \a» &y aW 

or, by the equation of continuity, 



1 JDlo 9 1 JD^O A 1-^ /rn __2\ a D / ">\ A 

^T57~^W =0 ' or f^ 1 '"*' ,)=0 ' or f/^" ,)=b0 ' 

This is the law of adiabatic expansion of a monatomic gas. The neglected 
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terme hence give the necessary corrections to this law, which are due to the 
causes indicated. The left-hand side represents the net rate of increase of 
energy of molecular agitation, being the increase corresponding to the rising 
temperature after allowing for the change due to adiabatic alteration of 
volume. This net rate is analysed as follows, taking the various terms 
in order : — 

(a) Work is done on the molecules by the external forces, in the motion of 
interdiffusion relative to the resultant stream velocity c . 

(b) The variation in the proportions of the molecules mi and m 2 , due to 
diffusion, produces a flow of energy of stream-motion c > since \m\<$ exceeds 
! m 2 c 2 . 

(c) The viscous forces cause dissipation of energy, represented by the 
third term on the right-hand side. 

(d) The fourth term represents the effect of conduction, together with a 
hitherto unrecognised term due to diffusion. Supposing, for simplicity, that 
there is no mass motion, the equation of energy becomes 



at J J ax J 



CO 



7 w 



8T f -ao 8 ' °° 

4- |-Bo-~ < 2(ny8r + ^2/8- r )— ^-~%r(via r -\- i/ 2 a- r ) 
o& L 1 «o 1 



on substituting the values of U1C1 2 and U2C2 2 , and substituting for f ' in terms 
of ?V and dTo/di\ If, further, there is no diffusion, it reduces to the equation 
of conduction of heat with thermal conductivity #, viz., 



p olo v ^ / o 0I0 



where 



BT> f ^ /3 ' * 1 

3 =e 'J-Bo-y^-S S(/']/3 r + ^ 2 /3-r) ^r^i^-,) f. 

•J Li «o 1 J 



It appears also that the motion of interdiffusion is accompanied by a flow 
of heat proportional to the velocity of diffusion, a process which we shall 
term the " thermal flux of diffusion/' In the absence of conduction, external 
forces, and mass motion, the equation of energy is thus seen to be 



ox ax 



where 3© is defined by 



CO 



PT 1 

—^ , Z r {V] a r -f- i>2&~ r ) i 

J «o 1, 



and is called the " specific energy of diffusion/' 
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(4) The Inequality of the Partial Temperatures.. 

In I (3) we found that 

To'^DoSo'ToSXo'M 
so that the temperatures of the component gases are unequal by an amount 
proportional to the rate of change of the ratio of their densities. We will 
write the equation in the form 

To' = -IA7&, Ia = -iDoSbTo, 

I x being termed the anisothermal constant of diffusion. Since I A is found to 
be positive, the gas of diminishing relative density is the hotter. 

Part III. — Discussion of the Eesults in Various Special Cases. 

(1) Introductory Remarks and Notation. 

The various results arrived at in Part II are complete and perfectly 
general. To be of service to natural philosophy, however, we must consider 
itheir values for selected particular molecular models, and in order to obtain 
numerical results we, must be prepared to make successive approximations to 
the accurate formulae. Our usual procedure in this Abstract will be to state 
the general form of the first approximation, and to determine the correction 
«due to further approximations only for special types of molecule. The types 
which will be dealt with are (a) point centres of force varying inversely 
«as the ^th power of the mutual distance; (6) rigid elastic spheres; and 
i(c) rigid elastic spheres surrounded by a field of attractive force. 

The following notation must first be explained, in order that the formulae 
of this section may be understood. 

fii = r/i 1 /(m 1 + 7n 2 ), fi 2 = m 2 /(mi-f m 2 ), fi\2 = mi/m 2 = /t*i//* 2 = I//&21, 



00 



<t>i2 k (y) = (2k 4- 1) {fiifiJh (% + m 2 ) } *jm { 1 — P* (cos 12 ) } p dp. 

Here 6\ 2 is the angle through which the relative velocity of two molecules, 
.mi, m 2 , is changed, during an encounter in which p is the distance 
between their initial or final lines of undisturbed rectilinear motion, and 
|(mi 4- m 2 )/ A mim 2 )] ^ 2/ is the magnitude of their initial or final relative 
velocity ; P# (cos i2 ) is the Legendre function of cos 6\ 2i of degree k. If 
both molecules are of the same kind, ?% or m% the function corresponding 
to $\ 2 k (y) is written <f>n k (y) or <f> 22 k (y). In the next equation, 

K * (0 = 771TTTS P"* 8 ** &> ij2it * k+1) (ly > 

\i-t-h-rji)t+kJo 
X /c (t) is affected by the same suffixes as the function <f> k (y) under the 
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integral sign ; K and <£ clearly depend upon the molecular model adopted,. 

the various models differing in respect of the dependence of Q\ 2 on p and y.. 

Again, 

h = Kia 1 (O/K12 1 (0), k® = K 2 (0/Kia' (0) ; 

in the latter the suffixes 11, 12, or 22, are to be added to correspond to the 
suffixes of K 2 (t) on the left. 

(2) Particular Molecular Models. 

If the molecules are rigid elastic spheres of radii 01, c 2 , it may readily be* 
shown that 

\ho7r7n>im 2 / 

Jc t = <t + 2)tl(t + %) u kn® = [2oi/(<ri + <r a )] 2 *ia', h 2 ^ = [2<7 2 /(<7 1 + <7 2 )] 2 ft 1 /, 

ftia<*> = |.(< + 3)* +1 /(^ + |) m . 

In these formulae (£ + 2)$, or, more generally, r«, where £ but not necessarily 
r is a positive integer, denotes r(r— l)...(r~£ + l). 

When the molecules are wth power centres of force, the force at unit 
distance being Ki 2 Wim 2 , it may be shown that 

Krf(O) = 2 I, W (lo»\ 2 ^ 1 >-* [Ki2(mi+m2)F) - J)r / 3 _^_ 
7r 3/2 / \m 1 + m 2 / \ n — . 



fo 



r p+3-2/(n-i) ] 



(« + l)»r [(3-2/(n-l)]' 
, w = 5 r[* + 4-2/fo+l)] I,(n) 

(^+D^ir[3~2/(^~i)]i 1 (72) ? 

ft n (0 = (Kn/Kia) 2 ^- 1 )^', W> = (Kaa/Kw^- 1 )^', 
where Ii(ti), Ia(^), are pure numbers, depending on % only.* 

When 71 = 5, so that the molecules are Maxwellian, we may note that 
k t = 1 for all values of £. 

When the molecules are rigid elastic spheres surrounded by a field of 
attractive force, the latter may be allowed for very approximately by an 
additional factor (1-fSia/T) in Ki^O), a factor 

{l + [2/0 + 2)]S ia /T}/(l + Sia/T) 

inland a factor {l + [3/(* + 3)]S/T}/(l + Si 2 /T) in U*\ Here S 12 , which is 
known as Sutherland's constant,! is given by 

Z>ia 2 mim 2 /4 (mi + m 2 ) K, 

* They are the same as the quantities thus denoted in •- Jeans' * Dynamical Theory of 
Gases/ 2nd ed., §§305, et seq. 

t S 12 is Sutherland's diffusion constant ; when the gas is simple, S 12 becomes S n , but 
Sutherland's constant of viscosity (generally written C) is f S n , or b^mJlZH. 
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where J&X2 2 is the potential of the force between two molecules m h m 2 , when 
in contact ; Sn and S22 are similarly defined, and appear in the numerators of 
the additional factors in Jen® and J^^- 

By putting m x = m 2 , "v% = 0, and changing the suffix 2 into 1 throughout, 
our formulae for a composite gas reduce to those appropriate to a simple gas. 
We shall first deal with the latter. 

(3) The Viscosity of a Simple Gas. 

The general first approximation to the coefficient of viscosity ku of a 

simple gas is 

5KT/™ii 2 (0). 

This result was oiven and discussed in an earlier memoir,* so that onlv the 
correction introduced by further approximations need be considered here. 

In the case of rigid elastic spherical molecules, the correction factor is purely 
numerical, as also in the case of n\X\ power centres of force. Successive 
approximations to the factor in the former instance are found to be 1*01485, 
1*01588, 1*01607, so that 1*016 times the above result may be taken as exact- 
to well within | per cent. For ?ith power centres of force, similarly, the 
correction factor is found to have the values 1*000, 1*004, 1*007, 1*011, 1*016, 
corresponding to the typical values 5, 9, 15, 25 and infinity for n. 

In the case of rigid elastic attracting molecules, the correction factor varies 
very slightly with temperature. Typical values, corresponding to the values 

00 , 5, 3, 1, 0*4, 0*2, 0*1, for S/T, are respectively 1*016, 1*010, 1*006, 1*000, 
1*004, 1*014, and 1*015. Thus, in all cases, the first approximation is 
exceedingly close to the exact value. 

(4) The Thermal Conductivity $n of a Simple Gas. 
A first approximation leads to the result 

f^ 

•as given in the earlier memoir just cited. As in the case of viscosity, 
the correction consists of a numerical factor to be applied to the above 
result ; this factor is very nearly unity, being 1*010 for rigid elastic spheres, 
and for nXh power centres of force varying from 1*000 for n = 5, through the 
values 1*003, 1*006, 1*007 for n = 9, 15, and 25, to 1*010 for n = go . 
Similarly for attracting spheres, as S/T takes the values 00 , 5, 3, 1, 0*4, 0*2, 

01, 0, the factor is found to equal 1*022, 1*007, 1*003, 1*000, 1*002, 1*003, 
1*004, and 1*010 respectively. The corrections to the first approximation are 

* 'Phil Trans.,' A, vol. 211, pp. 433-483, §§17, 20 (1911). 
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hardly within the present experimental powers of measurement, and the 
simple formula ^tcnO v agrees well with the observed values for argon, helium, 
-and neon. 

(5) The Coefficient of Diffusion Di 2 . 

The general first approximation to D12 is 

3 (mi -f m 2 ) ET/2 tt^o^i^K^ 1 (0). 

This result was first obtained by Langevin, and subsequently, independently, 
in the earlier memoir by the present author, already cited. Like the first 
approximations also to /en, 5n, and #12, it is an exact formula in the case of 
Maxwelliari molecules, and in no other. Maxwell obtained this expression 
for D12 in his own special case, though without a rigorous proof ; his method 
when applied to other cases (as by the authors named) leads only to the above 
first approximation, the error of which is in some cases considerable. We 
shall write the correction factor introduced on a second approximation in the 
form 1/(1 — eo). Before considering the value of eo we shall examine some 
particular cases of our formulae which throw light on the accuracy of a second 
approximation in general. 

When the mass and size of the molecules m 2 are negligible compared with 
those of mi, an exact solution for D12 is obtainable, as Lorentz first showed. 
The first few approximations to the correction factor in this case, when the 
molecules are rigid elastic spheres, are found to be 1*083, 1*107, 1*117, 
-converging to the value 32/97T or 1*132. If the molecules are nth power 
centres of force, the correction factor is equal to 



9\ 11 — l/n — 1\ n — ll\ n— 1/ sin[27r/(w — 1)]' 

which, when n is 5, 9, 13, 17, go , has the respective values 1*000, 1*031, 1*056, 
1*072, 1*132. These results agree with Lorentz's theory, and show that in this 
case, at any rate, the first approximation to D12 is seriously in error. 

The next simplest case of a general character with which we can deal is 
that in which the molecules mi, m 2 are identical in dynamical properties 
{e.g., mass, size, force), as, for instance, in the diffusion of a gas into itself. 
Hence we shall in this case w r rite Dn for Di 2 . The correction factor is 
found to depend only on the nature of the molecule ; its values for wth 
power centres of force, when n equals 5, 9, 17, or co , are respectively 
1*000, 1*004, 1*008, or 1*015, and the last value also corresponds to rigid 
spherical molecules. These are only first approximations, but would not be 
altered by further approximations by more than 0*001 or 0*002. The 
♦complete formula for rigid spherical molecules is 

0-1520/v! (2<70 2 (Aimi) 1 ^. 
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The correction factor in this case had been previously calculated approxi- 
mately by Pidduck, who obtained 0*151 (given thus to three decimal places 
only) in place of 0*1520 as above. 

For the same type of molecule, we also have Du = 1'200/cnfpi, which 
may be compared with an approximate result obtained by Jeans, using the 
" mean-free-path " method, viz., Dn = l*34tfn//oi. For other types of 
molecule, Dn is a different multiple of K U /pi. 

(6) Variation of Di 2 with the Ratio of Mixture v\\v%. 

The above first approximation to Di2, it may be noticed, depends only on 
v\ + v% and not at all on the ratio v\ : v% Meyer disputed the independence 
of D12 on the ratio of mixture, and developed a theory predicting considerable 
changes of D ]2 as v\\v% was varied. Experimental results showed a small 
variation of Di 2 , but much less than Meyer's formula alleged. This cast 
suspicion on both Maxwell's and Meyer's theories ; it now appears that the 
former was correct, but failed to represent the facts because of the unsuit- 
ability of the molecular model. A Maxwellian gas is the only one for which 
D12 is independent of v\\v 2 . The correction factor to the above approximate 
formula is in every other case a function of v\\v%. 

We shall consider only the approximation 1/(1 — e ) to the correction 
factor ; this is very nearly exact. We find that 

eo = (A/1 — i; -i — ytth 7~i — 2' 

d\v? + 2 di2Viv 2 + d 2 V2 
where 

&1 == 2 JULifJbvAl ', h = 2/^2^12^22°, 2& 12 = 30 (l — 4/Lt 1 /A 2 ) + 8^/^12°, 
fl2 

2d 12 = ~ — kn°k*° + i(l-4w)(l-^^ 

Zu fJb\fJb2 

and d 2 is obtainable from d\ by interchanging the suffixes 1 and 2. 

The behaviour of € with change of vijv 2 is easily studied when the 
molecules are rigid elastic spheres. In this case (remembering the 
convention mi^^) it appears that hjdi exceeds h 2 /d 2 , and that hi 2 jd X 2 
exceeds h 2 /d2. Further, if 

i(Tl(T2 "1 2 ^ 1 JU 2 3 



>i^l(l + 80 / i 2 ~100 / i 1 M 



.(o-i + cr 2 ) 2 J ^ 4 fn 

we have h/di ^> ^12/^12 ^> ^2/^2 ; the signs of equality correspond to the case 
mi = m 2 , ar\ = <T2, as in self diffusion. This indicates why Dn does not 
depend on vi/v 2 . In the general case, if the above condition is fulfilled, 



Simple and Composite Monatomic Gases. 



17 



€ and 1/(1 — €o) steadily decrease, as the proportion of molecules m 2 is 
increased, while, if it is not fulfilled, e and 1/(1 — eo) (and consequently also 
D 12 ) first increases and then decreases, as A 2 increases from to 1. Clearly, 
the more equal the molecular masses, the more equal also must he the 
molecular radii, in order that the above condition may be fulfilled. 

As regards the magnitude of the correction factor, it is greater the more 
unequal the diameters and masses of the molecules, reaching a maximum of 
32/97T, or 1*132, when mi/m 2 and 0-1/02 are infinite. This is for rigid 
spherical molecules : for other types the range is less, shrinking to zero in 
the case of Maxwellian molecules. In other cases the factor appears always 
to exceed unity. 

The magnitude of the variation in ordinary cases, and the comparison of 
experiment with theory, may be gauged from the following Table* ; the 
theoretical values of D i2 are so calculated as to make their mean agree with 
the mean of the experimental values, so that the theory is concerned only 
with the variations from the mean. The molecules are assumed to be rigid 
spheres, though a less simple model might prove still better. 



VTciSGS. 


A Vl 

At — , 


Correction 
factor. 


D12. 


Observed. 


Calculated. 

1 

i 

1 


Areon and helium 


0-273 
0-315 
0-377 
0-500 
0-677 
0-763 

0-25 

0-5 

0-75 


1-050 
1-056 
1-063 
1-073 
1-087 
1-094 

1-050 
1-073 

1-097 


0-244 
0-250 
0-250 
0*254 
0-256 
-263 

0*276 
0*280 
0*289 


i 
-248 
-250 ! 
-251 ! 
-254 ! 
-257 1 
-259 ! 

1 
j 

-276 ' 


Oxygen and hydrogen 




'282 j 
-289 

j 



It may be added that Kuenenf has improved Meyer's theory of D i2 by 
taking persistence of velocities into account, and has thus much reduced 
Meyer's predicted variations with v\\v 2y though they still remain twice or 
thrice as great as the observed changes. 



* The observed data were obtained in the Halle Laboratory, under Prof. Dora, by 

■n 

Schmidt, Deutseh, Jackmann, and Lonius (' Ann. d. Phys.,' vol. 29, p. 664, 1909). 
f Kuenen, Supp. 8 to 'Leiden Phys. Comm.,' January, 1913. 

VOL. XCIII. — A, C 
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(7) Thermal and Pressure Diffusion. 
The first approximation to the coefficient of thermal diffusion Dt is 

3(wii + m 2 )BT (k 1 -X)e 1 
27r7->o?%m 2 Ki 2 1 (0) z^oAi 
where 

6l = 30 (vifi 2 + v 2 m) (pi — /J*) — 8'fiifijci2° (yi —v 2 ) + 2 (yiftisJcn — ^2^21^22°), 
and 

Ai = {ljviv 2 ) { (cUvx 2 + 2d ia v 1 v 2 + d 2 v 2 2 ) — (h — l) 2 (hvi + 2b 12 viv 2 + hv£) }. 

In the case of Maxwellian molecules D T vanishes (as also every approximation 
to it, as here, since h\ = 1 for such a gas) ; like the variation of D 12 with 
v\\v% thermal diffusion is a phenomenon which wholly disappears in this 
alone out of all typical gases, so that the model is far from suitable 
physically, whatever its mathematical advantages. 

Thermal diffusion depends essentially on a difference of mass and structure 
between the molecules- mi, m 2 ; if m x = m 2 and Jen = k 22 °, etc., D T vanishes, 
as the above first approximation also suggests. The more unequal the 
molecules, the greater is D T . Considering elastic spherical molecules, if 
en = a 2 but m\jm 2 is very large, D T reduces approximately to 

D12 . 15X1X2/ (39 + X2) 

(thus if Xi = X 2 = J, D T = 0*095 D i2 nearly). When both mi/m 2 and ai/cr 2 
are very large (Lorentz s case) 

j\ _ 3X 2 Di 

where first, second, and third approximations to the numerical factor D//D' 
give 0*417, 0*494, 0*524, converging to 0*58 nearly — this also illustrates the 
degree of correction introduced on further approximations after the first, 
although the error of the latter is at its maximum for this specially unequal 
molecular pair ; the above formula is, moreover, true only so long as X 2 is not 
too nearly equal to 1, since D T vanishes with XiX 2 always. If Xi = X 2 = i, 
we get D T = 0*26 D i2 nearly, which may be compared with the former case 
when en = <r 2 , as illustrating the effect of the inequality between the molecular 
diameters when mi/m 2 is very large. 

Pressure diffusion likewise depends on the inequality of the molecules (of 
their masses only, however). Thus 

D^ = D i2 . XiX 2 (mi— m 2 )/(Ximi4"X 2 m 2 ). 

The following Table illustrates the relative magnitudes of D12, D T and D^, 
for a few typical gas-pairs (rigid elastic spherical molecules being assumed) : — 
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Katio of 
mixture, 


A-He 

m l /m 2 = 9 *8 o^/erg = 1 *7. 


O-H 

tni/ i m 2 = 16 ^i/^ = 1 *3. 


O-N 
mijm 2 =1-1 0*1/^2 — 1 '0. 


Dt / D 12 . 


IWDi* 


D^/Dt. 


Dt / D 12 . 


D//D 13 . 


D^/Dt. 


Dt/D 12 . 


r>/ / d 13 . 


IV Dt. 


1 : 3 
1 : 1 
3 :1 


0-133 
0-132 
0*079 


0-513 
0-408 
0*166 


3*9 
3-1 
2-1 


0-122 
0-128 
0-079 


0-592 
0-441 
0*230 


4-8 
3-4 
2-9 


0-008 

o-oio 

0-007 


0*026 
0-033 
0-024 


3-2 
3 3 
3-4 



In this connection we may also remark on the steady state of a gas 
without diffusion, under the influence of (a) external forces or (6) of a 
permanently maintained non-uniform temperature distribution. In the 
former case the pressure will, of course, be non-uniform. The equations of 
diffusion show that a permanently non-uniform composition will thus be set 
up, depending in case (a) on the ratio Di 2 7^i2 or D^/D^, and in (6) on 
D T /Di 2 . Both cases are exemplified in the upper strata of the atmosphere, 
above the region of convection, and in a future note the effect of temperature 
gradient on atmospheric constitution will be indicated in detail ; it is less, 
however, than the effect of the diminishing pressure. The influence of non- 
uniform temperature on composition has been verified experimentally, the 
proportions of two gases being found to differ in two bulbs kept at different 
temperatures and with an open connection between them ; a preliminary 
account of these experiments will appear shortly. The magnitude of the 
effect is easily indicated by an example ; if argon and helium, or oxygen 
and hydrogen, are the gases, mixed in originally equal proportions, and the 
hot bulb is kept at 100° 0. and the cold at 0° C, the heavier gas will con- 
centrate slightly in the cold bulb, 52 per cent, of it to 48 per cent, of the 
lighter gas, instead of the original 50 per cent, to 50 per cent., and vice versa 
at the hot bulb. 

(8) Viscosity and Conduction in Mixed Gases. 

The formulae for # 12 and 5 i2 in mixed gases are very complicated. A first 
approximation to #12 agrees with that given in my first memoir, but the first 
approximation to S12 was there given wrongly, since, in deriving it, it was 
assumed that conduction could take place in a uniform gas without diffusion, 
which we have just seen to be impossible. The correct formula is given in 
the second memoir here abstracted, but will not be reproduced in this place. 



(9) Inequality of the Temperatures of the Component Gases. 

This phenomenon has been shown to depend on dx® fit, or the time rate 
of variation of the relative proportions of the component gases. The following 
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example will illustrate its order of magnitude. Consider a gas composed 
initially of equal numbers of molecules which are similar to one another in 
mass and other dynamical properties (since the phenomenon does not, like 
thermal diffusion, depend on molecular differences, this assumption may 
legitimately be made, in order to obtain a simple numerical illustration). 
We will suppose that the gas is at normal temperature and pressure, and that 
the molecules are similar in mass and size to oxygen molecules, though 
differing from one another in some other respects. Then if the proportions 
are changing at such a rate that at the end of one second there will be 
51 per cent, of one component to 49 per cent, of the other, a temperature 
difference of about one thousand- millionth of a degree Centigrade will be set 
up, the hotter gas being that which is in diminishing ratio. Thus equi- 
partition of energy is disturbed only by a negligible amount even in this 
rather extreme case. 



On the Mechanical Relations of the Energy of Magnetisation. 

By G. H. Livens, The University, Sheffield. 

(Communicated by Sir J. Larmor, F.K.S. Eeceived September 19, 1916.) 

1. There still appears to be some uncertainty in the expression for the 
total energy of a magnetic field, and the mode of its separation into its 
fundamental constituents, as belonging to the aether and the matter. This is 
chiefly due to certain discrepancies of sign which exist in the results of the 
theory in its statical and dynamical aspects. In the statical theory the 
energy is usually made to appear as though distributed over the field with 
the density /aH 2 /87t, which corresponds to an sethereal density of amount 
H 2 /87r, but in the dynamical theory the same expressions are obtained with 
opposite signs. In some quarters it is considered that there are difficulties 
of a fundamental nature involved in any attempt to remove this dis- 
crepancy, and several authors* have tried to construct a more consistent 
theory on a new basis. 

The object of the present note is mainly to prove that it is possible to 
interpret the older and more usual form of the theory in a perfectly logical 
and consistent manner, so that the aforementioned discrepancy does not 

* Cf. ' Encyclopadie der mathematischen "Wissenschaften,' vol. 5, Art. 15, " Electro- 
statik u. Magnetostatik " (B. Gans), p. 338, where full references are given. 



